The structure of /o(980) from charmed mesons decays 
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Abstract 

We use the QCD sum rules to evaluate the form factors associated with the semileptonic decays 
of Dg and D mesons into /o(980). We consider the /o(980) meson as a quark-antiquark state 
with a mixture of strange and light components. The decay rates are evaluated in terms of the 
mixing angle. Using the same form factors to evaluate nonleptonic decays in the framework of the 
factorization approximation we conclude that the importance of the light quarks in /o(980) is not 
negligible. 
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Recent experimental data coming from charmed mesons decays have opened new possibil- 
ities to understand the spectroscopy of the controversial light scalar mesons, since they are 
abundantly produced in these decays 0. Actually, with the firm identification of the scalars 
a and k, the observed scalar states below 1.5 GeV are too numerous to be accommodated in 
a single qq multiplet. This proliferation of the scalar mesons is consistent with two nonets, 
one below 1 GeV region and another one near 1.5 GeV y. In this new scheme the light 
scalars (the isoscalars (t(500), /o(980), the isodoublet k and the isovector ao(980)) would 
form an SU(3) flavor nonet. In the naive quark model the flavor structure of these scalars 
would generically be: 
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= —j={uu — d(i) , = ag = du, (2) 



= US, = ds, = sd, K = su, (3) 

where we have already allowed a mixing between the isoscalars ss and {uu + dd) . Although 
the predominant ss nature of the /o(980) has been supported by the radiative decay cf) — *• 
/o (980)7, aiid by some theoretical calculations , there is no fundamental theoretical 
reason to expect a = 0, as we find in the vector meson sector. Actually, since instantons are 
supposed to be as important in the scalar sector as they are in the pseudoscalar sector j^, 
we would expect a mixing in the scalar sector similar to what we have in the pseudoscalar 
sector. In this sense, the measurements of J/ip /o(98O)0 and J/ip ^ /o(980)co' with 
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similar branching ratios |6|, indicating that /o(980) is not purely an ss state, can not be 
taken as a surprise, were used to estimate the mixing angle in 

Eq- ©, giving a = (34 ±6)°. A similar mixing angle, 35° < a < 55° j^, was found analysing 
the experimental results — > /o(980)7r"'' and — > ^vr"*" |lJ |. Using /o(980) as a pure ss 
state {i.e., using a = 0) the authors of ref. (ll|, could not reproduce the experimental result 
of Dj^ — > /o(980)7r+/ 07r+ [l^. They concluded that there is room for a sizable light 

quark component in /o(980), corresponding to a mixing angle of about a ~ 40°. 

In this work we propose that experimental results of the semileptonic decays — »• 
/o(980)£"'"z/£ and — > fo{980)i~^h'e, can be used to get the minimum bias estimate of 
the importance of the light quark content in the scalar-isoscalar /o(980), in the quark- 
antiquark scenario. The hadronic part of the current of the decay, can produce the 
/o(980) only through the ss component, while the hadronic part of the decay can produce 



the /o(980) only through dd. To observe these decays it would be necessary a high statistics 
experiment and also a tagging to separate the semi-leptonic from the decays. These 
two conditions could be satisfied next year by the CLEO-C experiment 

In order to estimate theoretically these semileptonic Dj^ and decays, we use the 
method of QCD sum rules which has been successfully applied to several semileptonic 
decay processes. Since the semileptonic decay is supposed to occur on the quark level, the 
decay rate should depend crucially on the direct coupling of the resonances to the quark 
currents. The coupling of the /o(980) to the scalar current 

js = ss cos{a) + [uu + dd) — ^— , (4j 

can be parametrized as 

{0\js\fo (p) ) = A/o = nif, ff^ , (5) 
and can be determined by the QCD sum rule based on the two-point correlation function 

n(g2) = ^ I d'xe^^-{0\T[j,{x)jmm. (6) 

This same correlation function was studied in ref. in the case a = 0. Following the same 
procedure and considering a general mixing angle a we get the sum rule (up to order m^): 
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In the numerical analysis of the sum rules, the values used for the strange quark mass 
and condensates are: = 0.14 GeV, {ss) = 0.8{qq), {qq) = —(0.23)^ GeV'^, {sgsa.Gs) = 
ml{ss) with ml = 0.8 GeV^. We evaluate the 2-point sum rules in the same stability window 
found in 14|: 1.2 < < 2.0 GeV^. In Fig. 1 we show the different contributions to Xj^ 
as a function of the Borel mass using the continuum threshold uq = 1.6 GeV^ and the 
mixing angle a = 37^. We see that A^^^ is very stable, as a function of the Borel mass, in 
the considered Borel range. The coupling is not very sensitive to changes in the values of 
the continuum threshold and mixing angle. For uq = (1.6 ± 0.1) GeV^ and < a < 37° we 
obtain 

A/o = (0.19 ±0.02) GeVl (8) 
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FIG. 1: Various contributions to the OPE of the couphng A^^^, as a function of the Borel parameter 
M^. The sohd Une gives the final result and the long-dashed, dashed, dot-dashed and dotted lines 
give the contributions of the first, second, third and fourth terms in Eq. ((T)). 

We can proceed and evaluate the /o(980) mass from the above sum rule, taking the derivative 
of Eq. (|7j) with respect to and dividing the resulting sum rule by Eq. ((7j). For uq = 
1.6 GeV^ and the mixing angle a = 37°, the mass is in a very good agreement with the 
experimental number and is very stable, as a function of the Borel mass. In the considered 
Borel range we obtain: rrif^ = (0.98 ± 0.01) GeV. For smaller (bigger) values of a we get 
a bigger (smaller) mass. The best agreement with the experimental result is obtained for 
a ~ 37°. 

The coupling given in Eq. (jHJ will be used in the analysis of the semileptonic decays of 
Ds and D into /o(980). The Dj /o(980)£z/^ form factors are defined through the matrix 
elements 

{W)\si,il - l,)c\DM) = cos(a) [f^^(t)ip + p'), + f^^{t)q,) , (9) 

and 

{/o(pOM7.(1 - l.)c\D{p)) = ^ {f^mp + p'), + fE{t)q,) , (10) 

with t = and q = p — p'. Since in the decay rate the form factor /-(t) is multiplied by 
the difference of the lepton masses, its contribution is negligible for both e and fi decays. 
Therefore, the differential semileptonic decay rates are given by 
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Xnml,ml,t){f^^{t))\ (11) 



for fo{980)i+Ui, and 

dT{t) sin2(a) Gl\Vj 



dt 2 1927r3mf) 



X'l\ml,mlMfm)\ (12) 



for D+ /o(980)£+^/^ In Eqs. (HH) and ^ A(x, ?/, ;z) = + + ^2 _ _ 2x2 - 2?/^, 
Gi? is the Fermi couphng constant and Vcs and are the Cabibbo-Kobayashi-Mashawa 
transition elements. 

Using the QCD sum rule technique the form factors in Eqs. © and ()1U|) can be 
evaluated from the time ordered product of the interpolating fields for Di and /o, and the 
weak current = qi-^^il — 75)0 (where qi is s or (i for Dj being Dg or respectively): 

7;(p,p') =i' j d'xdS (0|T[j,(x)j7(i/)jL^(0)]|0)e^(^'-+''-^) , (13) 

where the meson in the initial state is interpolated by the pseudoscalar current 

JdAx) = qi{x)i-i^c{x) , (14) 

and the /o(980) is interpolated by the scalar current given in Eq. Q. 

In order to evaluate the phenomenological side we insert intermediate states for Dj and 
/o, we use the definitions in Eqs. (jH)) and and obtain the following relations 

+ contributions of higher resonances , (15) 

for Dt /o(980) and 

rj^phen, _ ^Ifp . sin(a) /f'(t)(p + + f^{t)q^, 

+ contributions of higher resonances , (16) 

for D+ /o(980). 

In the above equations we have used the standard definition of the couplings of Di with 
the corresponding current: 

mn,m = ^^^l^. in) 

The three-point function Eq.lfT^ can be evaluated by perturbative QCD if the external 
momenta are in the deep Euclidean region 

p < (mc + rusf, p'^ < -iml, t < (m^ + m^)^ . (18) 



In order to approach the not-so-deep-Euchdean region and to get more information on the 
nearest physical singularities, nonperturbative power corrections are added to the pertur- 
bative contribution. In practice, only the first few condensates contribute significantly, the 
most important ones being the 3-dimension, quark condensate, and the 5-dimension, mixed 
(quark-gluon) condensate. For the invariant structure, {p + p')fj., we can write 



T 



+ T^='{q-iqi) + T^='{q-igsCT.Gqi) + ■■■ . 



(19) 



The perturbative contribution is contained in the double discontinuity p+. 

In order to suppress the condensates of higher dimension and at the same time reduce 
the influence of higher resonances, the series in Eq. p9|) is double Borel improved. Further- 
more, we make the usual assumption that the contributions of higher resonances are well 
approximated by the perturbative expression with appropriate continuum thresholds Sq and 
Uq. By equating the Borel transforms of the phenomenological expression for the [p + p')^ 
invariant structure in Eqs. (jl5p and and that of the "theoretical expression", Eq. (|19|). 
we obtain the sum rule for the form factor f^'{t). The sum rule for f+^{t) (at the order 
rris) is given by: 

«-2 -m2 _ ~^ 
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where C 



rric+ms 



^fo and 



p+{s,u,t) 



with sm 



A3/2(s,M,t) 



|n 2mcms{2ml — s — t + u) + ml{s — t + u) + s{—s + t + u) 
t + u){su + mcms{s — t + u))^ Q{s — sm), (21) 



m: 



+ ^2^^ - The sum rule for f^{t) can be obtained from Eq. ()20|) by just 
neglecting the rris terms, changing (ss) by (qq) and Dg by D. 

The decay constant foj defined in Eq. (fTTj) . and appearing in the constant C\ can also 



be determined by sum rules obtained from the appropriate two-point functions 
The value of the charm quark and meson masses are: rUc = 1.3 GeV, rufg - 
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0.98 GeV, 

1.97 GeV and mo = 1-87 GeV. For the Ds and D decay constants we use fo^ = 
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FIG. 2: Various contributions to the OPE of the form factor /^"(t), at zero momentum transfer, as 
a function of the Borel parameter M^. Sohd curve: total contribution; long-dashed: perturbative; 
dashed: quark condensate; dot-dashed mixed condensate contribution. 

(0.22 ± 0.02) GeV and fo = (0.17 ± 0.02] Gevjnl 15, Q- For the continuum thresholds 
we take the values discussed in refs. 11, [3, ll^h sq = (7.7 ± 1.1) GeV^ for Dg, sq = 
(6.0 ± 0.2) GeV^ for D and Mq = (1-6 ± 0.1) GeV. 

We evaluate our sum rules in the range 4.0 <M^< 8.0 GeV^ at a fixed ratio M^/M'^ = 
(m|,^ — m^)/ml , which is compatible with the Borel ranges used for the two-point functions 



in refs. |llL Il4| . In Fig. 2 we show the different contributions to the form factor f'^" at 
zero momentum transfer, from the sum rule in Eq. ((20)), as a function of the Borel variable 
M^. We see that the perturbative contribution is the largest one, and that the mixed 
condensate contribution is negligible. A similar behavior is also obtained for f^{t). Varying 
the continuum thresholds sq and uq, and the couplings A/^, fo^ and fo in the ranges given 
above, we get for the form factors at t = 0: 



(22) 



0.50 ±0.13 for D, 
0.53 ±0.15 for D 

Our values for f+'iO) are compatible with the values found in refs. 0, 0|, where the non- 
leptonic decay of the meson were studied. 

The t dependence of the form factors evaluated at = 7GeV^ in the range —0.5 <t< 
GeV^ can be fitted by a linear expression 



(t)=/f^^(0)±At. 



(23) 



This expression is consistent, in the considered range of momentum transfer, with a monopole 



expression f+'{t) = ^+ t , with the mass of the pole Mp = (1.70 ± 0.05) GeV for Ds and 



Mp = (1.95 ± 0.05) GeV for D. It is interesting to notice that Mp is compatible with the 

n n 

values found for the decays Ds f] [19|| and D ^ k jl8|. 

In the limits of the variables and the continuum thresholds discussed above we obtain for 
the semileptonic decay widths 

T{Dj /o(980)£+z/^) = cos2(a)(8.1 ±4.1) GeV , 
T{D+ /o(980)£+z/^) = sin2(a)(1.5 ± 0.8) GeV , (24) 

where we have used Vcs = 0.975 and Vcd = 0.22. The ratio between the two decay widths 
given in Eq. gives us a direct information about the mixing angle: 

r{Dt ^ /o(980)f+z/,) ^ 56 ±2 

r(D+ fo{980)i+ue) tan2(a) ' ^ 

For a ~ 37°, from Eq. ()25|) we can conclude that the semileptonic decay width of Dg into 
/o(980) would be around one hundred times larger than the semileptonic decay width of 
D into /o(980). Since in the semileptonic decays there are no complications due to strong 
interactions, we believe that an experimental measurement of the ratio in Eq. is the 
cleanest way to evaluate the mixing angle a. 

It is interesting to notice that if the current in Eq. (@]) were used in ref. for /o(980), 
instead of a pure ss current, the ratio calculated there would change to 

^- %:^rir ---w(o-^^^°-^B). (26) 

Therefore, for a ~ 37°, which is compatible with the mixing angle found in refs. 0,0], the 
above ratio would be reduced to i? = (0.28 ±0.11) in agreement with the experimental result 
j^exp ^ (0.210 ±0.069) 3. 

At this point we might conclude by saying that /o(980) is well described by the relation 
in Eq. (0), with the mixing angle being a ~ 35°, and by making the prediction in Eq. ()25j) . 
However, there are two experimental facts which do not fit in this picture and that might 
even be in contradiction with Eq.(fH): 

• A mixing angle in /o(980) also implies that cr would have a strange component (see 
Eq-O- Therefore, the decay Df an^ should also occur. As a matter of fact. 
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the ratio T{Df (T7r"'")/r(D+ /o(980)7r"'") is proportional to tan^(Q;). However, 
the E791 Collaboration did not observe any contribution from cxtt^ to the 

Df — > TT+TT+TT^ decay [3] pointing towards a ~ 0. 

In the framework of generalized factorization, the amplitude for the /o(980)7r"'" 
decay is given by (neglecting the annihilation term) 



A{D^ /o(980)7r+) = U{f^) + ^) - mjf^f^iO) , (27) 



since (0|V^|/o) = due to charge conjugation invariance and conservation of vector 
current. In Eq. (j27j) Cj is the Wilson coefficient entering the effective weak Hamiltonian. 
Therefore, using the results in Eq. we get 

r(Z}+ ^ /o(980)7r+) 46 ±2 



T{D+ -> /o(980)7r+) tan\a) ' 

From the E791 Collaboration, the experimental value of this quantity is jlol |: 

fr{D+ /o(980)7r+)y'^^ 



(28) 



Vr(D+ ^ /o(980)7rH 
which would lead to a ~ 62°. 



13.3 ±0.4, (29) 



It is important to remember that in the analysis of the nonleptonic decays we are using 
the factorization approximation. Therefore, the apparent inconsistency between the data 
and the mixing angle can still be due to this approximation. This is why the measurement 
of the ratio in Eq. would bring important information about this puzzle. If the light 
scalars could not be seen in the semileptonic decays, this would clearly indicate that their 
structure is more complicated that simple quark-antiquark states. One possibility is that 
they are four-quark states, as suggested in refs. j^l^^Isil- 

To summarize, we have presented a QCD sum rule study of the and D semileptonic 
decays to /o(980), considered as a mixture of the scalars ss and uu + dd. We have evaluated 
the t dependence of the form factors f+'{t) in the region —0.5 < t < OGeV^. The t 
dependence of the form factors could be fitted by a linear form compatible, in the studied 
range, with a monopole form, and extrapolated to the full kinematical region. 

The form factors were used to evaluate the decay widths of the decays /o(980)£+z/£ 
and D'^ /o(980)£"'"z/^ as a function of the mixing angle. Experimental data about these 
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decays would provide a direct estimate of the mixing angle. Using the same form factors 
to evaluate nonleptonic decays in the framework of the factorization approximation, it was 
not possible to explain all available experimental data with a fixed mixing angle. However, 
data seem to suggest that there is a sizable non strange component in the /o(980) meson. 
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